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I. INTRODUCTION 

Dynamical chiral symmetry breaking is certainly one of the most interesting phenomena 
of QCD. It is entirely a strong coupling effect in the sense that dynamical quark masses 
cannot be generated at any order in perturbation theory. Thus nonperturbative methods 
like lattice Monte-Carlo simulations P|, chiral perturbation theory or the Green’s 
function approach using the Schwinger-Dyson and Bethe-Salpeter equations (SDE/BSE) 
i9 , are needed to explore chiral symmetry and its breaking pattern. 

With the rediscovery of the Ginsparg-Wilson relation and the construction of actions 
for overlap, domain wall and perfect fermions, the lattice formulation of QCD emerged 
in principle as an appropriate nonperturbative tool, with which to study the effects of 
dynamical chiral symmetry breaking. However in practice, lattice simulations with small 
quark masses are extremely expensive in terms of CPU-time. It is only with staggered 
fermion actions that quark masses not far from their physical values have been achieved to 
date, but these actions have the disadvantage that full chiral symmetry is only recovered 
in the continuum limit. Consequently, there is no certainty with any hnite volume that 
the correct breaking pattern can be observed. 

Lattice simulations are, of course, always performed at a hnite volume. Since continuous 
symmetries cannot be spontaneously broken at a hnite volume V, chiral symmetry is 
restored in the limit of zero quark mass, m —0, independently of the formulation of the 
lattice action. Thus one hrst has to perform the limit U —>■ oo before one can investigate 
the chiral limit. At the hadron level, whether for mesons and baryons, chiral perturbation 
theory provides a reliable tool with which to make such an extrapolation (see e.g. 9901 
and references therein). Volume ehects for any particles that couple to the pion can be 
arranged in powers of exp[—M^rL], where M^- is the pion mass and L is the size of the box 
9]- On the other hand, chiral perturbation theory has nothing to say about volume effects 
in the underlying quark and gluon substructure. For this the Green’s function approach 
employing Schwinger-Dyson equations provides a suitable alternative. Indeed, a recent 
investigation of the gluon and ghost propagator of the Landau gauge Yang-Mills theory 
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on a torus could offer an explanation for systematic differences in the infrared behaviour 
of these propagators in a box, compared with the inhnite volume limit . 

In the present work we extend this analysis and investigate hnite volume effects in the 
quark propagator. To this end we study the QCD-gap equation for the quark propagator 
on a torus. The input from the Yang-Mills sector of QCD consists of numerical solutions 
for the ghost and gluon propagator, which match corresponding lattice calculations ET 


1 1| |. Presently unknown contributions from the quark-gluon vertex are parametrised such 


that the quenched lattice quark propagators from each of Ref. [l^ (staggered) and Ref. m 
(overlap) are reproduced by the gap equation. This idea has already been explored to 
some extent by Bhagwat et al. (l4l |. 

Our treatment differs from that of Ref. in two essential respects: hrstly, we solve 


the coupled set of three SDEs for the ghost, gluon and quark propagators, and so include 
the Yang-Mills sector of the SDEs in the Green’s function approach. Secondly, we calcu¬ 
late our propagators on a manifold similar to that of the lattice and £t the interaction 
to reproduce the lattice results of the gluon and quark propagators on their respective 
manifolds. 

One of the advantages of the Green’s function approach is that volume effects can be 
studied continuously from very small to very large volumes (corresponding studies for 
meson observables using chiral perturbation theory e.g. have to distinguish between two 
different regions of chiral counting 0 )- Furthermore one has direct access to the inhnite 
volume and the continuum limit without the need to perform any extrapolations. We are 
thus in a position to study chiral symmetry restoration at small volumes together with 
effects at large and inhnite volumes in the same framework. 

The paper is organised as follows: In the next section we shortly review basic properties 
of the pattern of dynamical chiral symmetry breaking in a box. In particular we recall 
the derivation of the Gasher-Banks relation and a basic estimate for a minimal box length 
for chiral perturbation theory. In section IIIII we discuss the technical details associated 
with solving the quark Schwinger-Dyson equation on a torus. Our results on the compact 
manifold are presented in section IIVI We start with a summary of hnite volume ehects 
in the Yang-Mills sector based on the results of Ref. . These are used as input into the 
quark-SDE together with a model for the quark-gluon interaction, which is discussed in 
The interaction is fitted so that the lattice data for the quark propagator 
(staggered) and Ref. US (overlap) respectively are reproduced by the quark- 
SDE. We then determine the corresponding quark propagator at larger volumes and in the 
inhnite volume limit. The results are compared in section HVBl For small quark masses 
we hnd sizeable quantitative volume ehects which are still present at comparably large 
volumes. These ehects are qualitatively similar for both the propagators from staggered 
and overlap quarks. The small volume behaviour of the quark propagator is investigated 
in section llV Gl For a hxed, small current quark mass we determine the onset of dynamical 
chiral symmetry breaking when the volume of the box is increased. We hnd a minimal box 
length of L ~ 1.6 fm, below which chiral perturbation theory cannot be safely applied. 
Finally, we discuss the inhnite volume properties of our quark propagator. We determine 
the chiral condensate in section IV Al comment on possible analytic structures in section 
IV Bl and give results for the corresponding pion mass and decay constant in section IV Gl 
We summarise and conclude in section ED 



from Ref. 
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II. CHIRAL SYMMETRY BREAKING IN A BOX 


Before we embark on our investigation, let us recall the finite volume behaviour of the 
chiral condensate, as this is the order parameter of dynamical chiral symmetry break¬ 
ing, 0. The fermion propagator in its spectral representation is given by 


SA{.x,y) 



Un{x) ul{y) 

m — iX ’ 


( 1 ) 


where Un{x) and A„ are eigenfunctions and eigenvalues of the Euclidean Dirac operator, 
^n{x) = XnUn{x). The gauge field A is treated as an external field. These eigenfunctions 
occur either as zero modes or in pairs of opposite eigenvalues. Setting x = y, integrating 
over X and neglecting the zero mode contributions, one obtains 


1 
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m^ + X^' 


( 2 ) 


The quark condensate can be deduced by averaging the left hand side of this equation 
over all gauge held conhgurations and then taking the inhnite volume limit to give 


(gq) = -2m [ dX —(3) 
J m^ -|- A"' 

0 


where p(A) is the mean level density of the spectrum, which becomes dense in the inhnite 
volume limit. In the chiral limit, m —0, only the infrared part of the spectrum contributes 
and one hnally arrives at the Banks-Casher relation 0 

{qq) = -7rp(0) . (4) 

If the two limits are interchanged, i.e. if one takes the chiral limit before the inhnite 
volume limit, one has a discrete sum in Eq. Q and the infrared part of the spectrum 
cannot trigger a nonvanishing chiral condensate: chiral symmetry is restored. If, however, 
at a given volume the explicit quark mass m is not too small, one can still observe the 
spontaneous formation of a quark condensate. If the factor (m^-|-A^)“^ varies only slightly 
with n, the sum in Eq. m can still be replaced by an integral and Eq. © remains valid. 
For this to be a legitimate approximation one needs m 3> AA l/Vp{X)=7r/{V\{qq)\), 
at the lower end of the spectrum. Thus one obtains the condition 


Vm\{qq) \ S> vr. 


( 5 ) 


To get a feeling for this condition, note that for a typical value of the chiral condensate 
of I (?*?) I = (0.25 GeV)^ and a volume of E = (5 fm)"^ the quark mass has to be of the 
order m ^ 5 ■ 10“^ GeV, which is well satisfied for all quark masses of physical interest. 
Whether there are sizeable modifications to the corresponding quark propagator due to 
the box is however a different question, as we shall see in section IIV Bl 
Ghiral perturbation theory builds upon the chiral limit, i.e. it can only be applied on 
volumes large enough such that small quark masses remain accessible. Gorrespondingly 
the chiral expansion parameter has to be small. On a torus the bosonic degrees 

of freedom have momenta p = 27rn/L with n a vector of integers. Small nonzero momenta 
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are therefore only present if the condition 


L » 


1 

w. 


1 fm 


( 6 ) 


is satisfied. A priori there is no way to say by how much L has to exceed 1 fm [^. In 
section ITV (II we will see that this scale can be estimated using the quark-SDE on a torus. 


III. THE QUARK SCHWINGER-DYSON EQUATION ON A TORUS 


In Euclidean momentum space, the renormalised dressed ghost, gluon and quark prop¬ 
agators in the Landau gauge are given by 






S{p) 


[, Pl^Pu\ Z{p‘^) 

J p2 

Zfjp^) 

i^ -|- M{jP) 


(7) 

( 8 ) 

(9) 


Here the ghost dressing function GijP'), the gluon dressing function Z{p^) and the 
quark wave function renormalisation Zf{p^) also depend on the renormalisation point 
/i^, whereas the quark mass function M{p^) is a renormalisation group invariant. These 
propagators are given by their corresponding Schwinger-Dyson equations shown diagram- 
matically in Fig. ^ Since we aim to analyse a quenched lattice quark propagator we will 
also work in the quenched approximation in the SDE approach and so neglect the quark 
loop in the gluon-SDE. An estimate of unquenching effects for the propagators and for 
light meson observables can be found in Ref. 0 . 

On a compact manifold, the ghost, gluon and quark helds have to obey appropriate 
boundary conditions in the time direction. These have to be periodic for the gluon and 
ghost helds^, and antiperiodic for the quarks. It is convenient, though not necessary, 
to choose the same conditions in the spatial directions^. We choose the box to be of 
equal length in all directions, Li = L 2 = -^3 = L 4 = L, and denote the corresponding 
volume V = L^. Together with the boundary conditions this leads to discretised momenta 
in momentum space. Thus all momentum integrals appearing in the Schwinger-Dyson 
equations are replaced by sums over Matsubara modes. Since the ghost and gluon SDE 
on a torus have been investigated in detail in Refs. we only discuss the quark-SDE 

here. On the manifold the quark-SDE can be written as 

s-'(p) = f d‘kj,s{k)r,{k,p)D,,(p-k), (10) 

where the factor Cp = 4/3 stems from the colour trace and we have introduced a reduced 
quark-gluon vertex ry(/c,p), by dehning rl^’'\k,p) = igY^u{k,p). The bare quark prop¬ 
agator is given by [5'°(p)]“^ = ij ■ p + Z^m^fj,'^), where is the renormalised current 


^ The condition for the ghost field can be read off easily from its BR ST-tran sformation |3. 
^ Different conditions in the spatial directions have been explored in Ellli 
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FIG. 1: The coupled set of Schwinger-Dyson equations for the gluon, ghost and quark propagators. 


quark mass. The wave function and quark mass renormalisation factors, and are 
determined in the renormalisation process. The ghost renormalisation factor, Z 3 , will be 
discussed below, when we introduce our expression for the quark-gluon vertex. The quark 
mass function M{p^) and the wave function Zf{p^) can be extracted from Eq. (fUHl by 
suitable projections in Dirac-space. 

Note that the quark propagator determined from Eq. m is independent of the regular- 
isation procedure. In our numerical calculations we use a subtracted version of Eq. m 
and an 0(4)-invariant UV-cutoff (for details see e.g. ref. [13 )■ is a simple matter to 
explicitly verify numerically (and also analytically) that the resulting quark propagator 
is independent of the cutoff, which therefore can be sent to inhnity at the end of each 
calculation. The quark-SDE, Eq. m, therefore represents not only the inhnite volume 
limit but also the continuum limit (in coordinate space) of any representation of the SDE 
on a compact manifold. We will use the phrase inhnite volume/continuum limit in the 
following to indicate this simultaneous removal of both an ultraviolet and an infrared 
cutoff. 

On a torus with antiperiodic boundary conditions for the quark helds, the momentum 
integral changes into a sum of Matsubara modes. 


/ 




jj E (■■■). 


711,712,ns,714 


( 11 ) 


counting momenta qn = X)i=i 4(27i'/L)(nj-|-1/2)6*, where e* are Cartesian unit vectors in 
Euclidean momentum space. For the numerical treatment of the equations it is convenient 
to rearrange this summation so that it represents a spherical coordinate system j^, see 
Fig.H for an illustration. We then write 
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FIG. 2: Two-dimensional sketch of the momentum grid dual to the four-torus for a fixed Cartesian momentum 
cutoff. The hyperspheres depicted by dashed lines are not complete in the sense that additional momentum points 
on these spheres are generated if the cutoff is increased. The 0(4)-invariant cutoff used in our calculations sums 
only over complete hyperspheres, which are indicated by fully drawn circles. 


^ E (■■•) (12) 

ni,712,^3,724 

where j counts spheres with qnqn = const, and m numbers the grid points on a given 
sphere. The corresponding momentum vectors are denoted qmj and their absolute values 
are given by qmj = |qm,j|- It is then a simple matter to introduce an 0(4)-invariant 
cut-off by restricting j to an interval [1,N]. The resulting quark SDE is given by 




N 






Note that the momentum argument of the gluon propagator is a difference of two antiperi- 
odic Matsubara momenta and thus lives on a momentum grid corresponding to periodic 
boundary conditions as it should. 

The quark Schwinger-Dyson equations, Eqs. m and m, can be solved numerically 
employing well established methods once the input from the Yang-Mills sector, the gluon 
propagator and the fully dressed quark-gluon vertex T^{k,p) are specihed. Our 
numerical method on the torus is outlined in Ref. , the corresponding continuum method 
as well as details on the renormalisation procedure of the quark-SDE a re g iven in Ref. 

The truncation scheme of the Yang-Mills sector is discussed in Refs. (2fl| |. 
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IV. NUMERICAL RESULTS FOR THE PROPAGATORS 
A. Yang-Mills sector and parameter fitting 

Before we discuss our ansatz for the quark-gluon vertex, let us shortly summarise the 
results of Ref. for the ghost and gluon propagators on the torus. In Fig. 13 the numerical 
solutions for the gluon propagator Z{p^)/p^ (left diagram) and the ghost dressing function 
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FIG. 3: The results for gluon and ghost from Dyson-Schwinger equations in the continuum and on the torus are 
compared with the lattice data of refs. min 


G{p^) (right diagram) in the continuum and on a torus^ are displayed together with the 
results of recent lattice simulations. Overall there is very good agreement between the 
DSE-solutions on the compact manifold and the lattice data. However, the infrared (IR) 
behaviour of the DSE-solutions on is qualitatively different, although the truncation 
scheme used in solving the SDEs is the same. The ghost dressing function in Fig. |21 diverges 
in the inhnite volume/continuum limit, whereas it stays hnite on the compact manifold. 
For the gluon propagator, this difference can be expressed in terms of an infrared power 
law, 

Z{p^) ~ ( 14 ) 


One obtains k ~ 0.5 (IR-£nite) on a compact manifold (even for very large volumes'), 
whereas k 0.596 (IR-vanishing) on R in agreement with analytical results |20L 1221 l20L 
23. This is a decisive difference, since it can be shown that an infrared vanishing gluon 
propagator cannot have a positive dehnite spectral function and is therefore conhned. 
Indeed, Zwanziger has argued that the lattice gluon propagator should vanish in the 
continuum limit and therefore be conhned as an effect of the proximity of the Gribov- 
horizon for low momentum gauge held configurations. However, no statement could be 
made as to the rate with which the continuum limit behaviour is approached. Current 
extrapolations of lattice data to the inhnite volume limit (on large asymmetric lattices) 
are somewhat ambivalent. Whereas an extrapolation of the gluon dressing function leads 
to K ~ 0.52 (IR-vanishing), one obtains k = 0.5 (IR-hnite) from an extrapolation of the 
corresponding gluon propagator 13 Therefore it seems as if there is a genuine diherence 
between propagators on diherent manifolds, which has to be taken into account when 
extrapolating to inhnite volumes. The diherences shown in Fig. Elmay serve as a measure 
of the upper limit of these ehects. 

The (quenched) solution for the gluon propagator on both type of manifolds is used 


® There are two slight changes compared to the treatment in Ref. |^: we adapted the overall scale to match the lattice 
results and the ultraviolet cutoff has been increased from A = 2.47 GeV to A = 3.60 GeV in order to minimise artifacts 
due to the cut-off. 







directly as input in the corresponding quark-SDE on a torus and in the inhnite vol¬ 
ume/continuum limit. What remains then is to specify an explicit expression for the 
quark-gluon vertex. Here we follow the strategy of Ref. □ and employ a parametrisa- 
tion of the vertex such that the (quenched) lattice results for the quark propagator are 
reproduced by the gap equation. Our ansatz for the vertex is 


= 7.ri(fc2)r2(A;2,/i2)r3(A;2,/i2) 


(15) 


where /c^ is the gluon momentum and is the renormalisation scale. The ansatz depends 
on the gluon momentum only and is thus the simplest possible form that respects charge 
conjugation symmetry. Furthermore this choice of momentum dependence ensures the ex¬ 
istence of a corresponding kernel in the Bethe-Salpeter equation for mesons in accordance 
with the axial-vector Ward-Takahashi identity, cf. section IV 01 The three components of 


this ansatz are given by 

ri(fc^) = 

’^7m 

(16) 

ln(fc 7 A^^^ + r) ’ 

^2(A(^/i^) = 

G(fc^/r2) G(7,/i^) ^3(h^) h [ln(fcVAgr)]^+^ 

(17) 

T,{k\p^) = 

a{M) -I- fc^/Ag^^ 

(18) 

where 5 = —9/44 is the (quenched) one-loop anomalous dimension of the ghost, 7 ^ = 
12/33 the corresponding anomalous dimension of the quark and r = e — 1 acts as a 


convenient infrared cutoff for the logarithms. It is well known that the effective interaction 
g^ZTiT 2 T 2 , in the quark-SDE has to approach the running coupling in the ultraviolet 
momentum regime ji^. In our ansatz this UV-part of the interaction is represented 
by Ti. The scale Aqcd is scheme dependent. Here, since we £t to the lattice data 
of Refs. EQ, its value corresponds to the MOM-scheme used therein. The product 
ZT 2 T^ goes to a constant for large momenta, since the ultraviolet behaviour of the ghost 
and gluon dressing functions is given by 


G{z) = G{s) a;log^^j-hl , Z{z) = Z{s) a;log^^j-Fl 


(19) 


with uj = / 3 oQ;(s)/( 47 r) = IINcOi{s) /and a large scale s. The anomalous dimensions 
of the ghost and the gluon dressing functions in ZV 2 T 2 , combine to 7 -|- 5 = —1 — 5, 
which is balanced by the explicit logarithm in r 2 . The scale Ag in Eq. (El represents 
a possibly scheme dependent scale inherent in the Yang-Mills part of the interaction, 
which is related to the analytic structure of the gluon propagator. The coefficient h is 
hxed such that the ultraviolet behaviour of the resulting running coupling matches the one 

The renormalisation group invariant G(C ^5 
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calculated on the lattice in Ref. 
with the arbitrary scale C, is introduced to impose the correct cutoff- and renormalisation 
point dependences of the effective interaction in the quark-SDE. Together, the product 
rir 2 represents the non-Abelian content of the quark-gluon vertex as expressed in its 
Slavnov-Taylor identity (STI) given by IIT 


G {k )KVy{q,k) = S (p) H{q,p) - H{q,p) S (q), 


( 20 ) 


where q and p are the quark momenta. This identity enforces the presence of the ghost 
factor G^k"^) in r 2 , which makes the quark-gluon vertex an infrared singular object, sim¬ 
ilar to the three- and four-gluon vertices [ 23 . The remaining part of rir 2 is infrared 
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staggered 

overlap 

^ia££jce(GeV) 

0.028 

0.057 

0.114 

0.090 

0.140 

0.210 

0.300 

MsDE{C){GeV) 

0.044 

0.080 

0.151 

0.076 

0.112 

0.162 

0.225 


TABLE I: Renormalised quark masses in the SDE on the torus at ^ = 2.9 GeV compared with the bare quark 
masses on the lattice from Ref. Ellll. Note that the bare staggered quark masses are smaller than the renor¬ 
malised quark masses from the SDE, whereas the corresponding bare overlap quark masses are larger (c.f. the 
discussion in the text). 



h 

AAGeV) 

Aqcd(G eV) 

ai 

a2 

03 

staggered 

1.33 

1.50 

0.35 

25.30 

4.80 

-1.39 

overlap 

1.31 

1.50 

0.35 

25.58 

3.44 

2.23 


TABLE II: Parameters used in the vertex model, Eqs. I15I18II . 


finite and can be interpreted as a model for the ghost-qnark scattering kernel H{q,p). 
The dependency of the vertex on the qnark wave function Zf through the inverse quark 
propagators in the STI is taken care of by Ts, the form of which is chosen appropriately. 
The extra factor Z 2 is vital in ensuring multiplicative renormalisability of the quark-SDE^. 
The dependence of this part of the vertex on the quark mass is expressed in terms of the 
function 


a(M) 


Oi 

1 + 02M(C2)/Aocd + osM 2 (C")/A|cd' 


( 21 ) 


In order to preserve multiplicative renormalisability of the quark-SDE, it is important 
that the scale ( at which the quark mass function is read off (and also the ghost factor 
in Ti) is not correlated with the renormalisation point. Instead it should be a hxed 
scale sufficiently far into the ultraviolet region that volume effects are negligible. In our 
calculations we use ( = 2.9 GeV. 

To £t the various parameters in our model interaction we solve the quark-SDE on a 
torus em ploy ing a momentum range similar to that used in the lattice calculations of 
Ref. EliJ. In the SDE this corresponds to a 24^ (36^) lattice in momentum space 
with a smallest momentum of 304 MeV (200 MeV) corresponding to a box-size of Li = 
L 2 = T 3 = L 4 = L = 2.04 fm (L = 3.10 fm) for the staggered (overlap) quarks. We 
first determine the parameters h,Ag and Aqcd by htting the ’Yang-Mills part’ TiT? of 
our vertex model to the ultraviolet part of the lattice running coupling from Ref. (ll| . 
which agrees with the two-loop results from perturbation theory. This fit is not unique 
but gives a range of pairs (A^jAgcm) related to the value of h. We then choose current 
quark masses m{pi?) such that the ultraviolet behaviour of the lattice quark mass functions 


^ Note that the effective interaction of Ref. Hi fails to ensure multiplicative renormalisability and is therefore only valid 
at the fixed renormalisation point chosen in their work. 
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are reproduced by the solutions of the quark SDE. As the ultraviolet behaviour of the 
quark mass function is determined by the ultraviolet behaviour of the input interaction, 
we hud the same result for each pair [Kg^KqcD)- This merely reflects the fact that the 
ultraviolet behaviour of the quark mass function is controlled by resummed perturbation 
theory and is therefore model independent. We then deduce values of the function a{m) 
in Eq. m such that the lattice quark mass functions are reproduced. Although this is 
possible for all values of the pairs {Ag, Aqcd), the corresponding quark wave functions 
Zf{p^) favour a small range of values for Aqcd and Ag. Finally we £x the coefficients 
ni, 2,3 and determine the values of the quark mass function M{C,) corresponding to different 
renormalised current quark masses Our hnal best parameter sets together with 

the quark masses are tabulated in Tables Him Note that the bare staggered quark masses 
are smaller than the renormalised quark masses from the SDE, whereas the corresponding 
bare overlap quark masses are larger. Comparing the functions Ms de{^ lattice) for both 
formulations we hnd agreement if the bare staggered quark masses are multiplied by a 
factor 1.70. No additional additive corrections, which recently have been identihed as a 
consequence of taste symmetry breaking 0 , 


are necessary. 


A few comments on the quality of our hts are in order. The scales Aqcd and Ag cannot 
be determined very well; we estimate the error in these scales to be of the order of 30 %. 
With the given values for these scales, the uncertainty in the parameter h is then only 

i.e. of the order of 
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related to the (small) error bars of the lattice coupling from Ref. 

2 %. The errors for our values of the parameters Oi, 02 , 03 are correlated to both the error 
in the determination of the current quark mass from the lattice data and the error bars of 
the lattice data in the infrared momentum regime. This leads to uncertainties of the order 
of a few percent. We have also performed hts to a(M) setting 03 = 0, thereby testing the 
possible redundancy of this parameter. The hts become worse in the range where lattice 
data are available, but still might be of tolerable quality. However, there is a nontrivial 
ehect which convinced us of the importance of the parameter 03 : if 03 = 0 we obtained 
vastly diherent values for a{M) in the chiral limit for the staggered and the overlap data. 
However, fitting 03 as a free parameter these values almost coincide. As a consequence, one 
obtains almost identical quark propagators if the inhnite volume/continuum limit and then 
the chiral limit are performed. On general grounds the overlap and staggered formulation 
of quarks on the lattice should coincide in the continuum limit. Possible residual scheme 
dependences at hnite current quark masses necessarily vanish in the chiral limit, provided 
there are no additive mass corrections (see above). This is exactly what we see if 03 is 
included as a free parameter. All qualitative conclusions of our paper do not depend on 
the details of the hts and are stable with respect to a variation of the parameters within 
the errors given above. For our quantitative results we have included these uncertainties 
in our estimate of the overall error margins. Unfortunately it is extremely difficult to give 
a quantitative estimate of the systematic truncation error, i.e. an estimate of possible 
changes once sub leading structures of the vertex are taken into account. Wherever possible 
we have tried to assess such ehects on a qualitative basis. 


B. The quark propagator at finite and infinite volnme 

In Fig. m we compare our results from the SDE with the lattice data. The solutions 
for the quark mass function M(p^) from the SDE can be nicely htted to the lattice 
results in both, the staggered and the overlap formulation. For the quark wave function 
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FIG. 4: The results for quark mass function M and the wave function Zf from Dyson-Schwinger equations on the 
torus compared with lattice data for staggered quarks m (upper panel) and overlap quarks [H (lower panel). 
The lattice momentum p has been used for the lattice mass functions, whereas the choice of the ‘kinematical 
momentum’ q corrects for hypercubical artifacts in the lattice wave functions min The results for the Ball- 
Chin ansatz, Eq. (HU, are shown only for the largest and smallest masses of the overlap quark propagator. 


Zf{p^) in the staggered formulation, shown in the upper right diagram of Fig. HI we 
observe a slightly larger spread in the SDE solutions than on the lattice®. The same is 
true when compared to the overlap data, although here in addition we observe a larger 
fall in the infrared. We have tried to reproduce this steeper decrease by considering 
various modihcations of our model interaction, Eqs. (Hii-dlHD, but have not succeeded. In 
particular introducing additional ghost factors, which correspond to an even more singular 
vertex in the continuum, does not improve the situation. It seems to us, that such a fall 
can only be reproduced when further tensor structures of the quark-gluon vertex is taken 
into account. Indeed, employing the Ball-Chin construction [Slf 

= i (A{p^) + A(i)) + 1 (A(/) - Mi)) + 

+ GB(pb-B(9b) (22) 


^ The parametrisation used in Ref. Hi leads to a similar large spread (not explicitly shown in their paper). 
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FIG. 5: The results for quark mass function M and the wave function Zf from Dyson-Schwinger equations in the 
continuum compared with the same lattice data as in Fig. 4 (staggered quarks in the upper panel, overlap quarks 
in the lower panel). Also shown is the chiral limit (dashed curves). 


instead of ji, (and modified parameters 01 ^ 2 , 3 ), we could at least reproduce the sharper 
decrease at small masses. However, as can be seen in Fig. the problem remains for 
large masses. It is therefore not clear to us, whether the sharper low momentum decrease 
in the overlap data should be taken seriously and really interpreted as an indication of 
the importance of a richer tensor structure in the quark-gluon vertex, or whether one 
should prefer the staggered data and conclude that the 71 ,-part of the vertex is sufficient 
to reproduce the lattice results. We leave this question open for future investigations in 
both the SDE-formalism and on the lattice, and proceed using our simple construction, 
Eqs. ffT3j) - (fTH|) . 

To assess hnite volume effects we now compare the lattice/SDE-results on the compact 
manifold with the inhnite volume/continuum limit. To this end we solve the quark- 
SDE in the continuum employing our lattice inspired ansatz for the quark-gluon vertex 
and the continuum solutions for the ghost and gluon propagators, discussed in section 
WM as input. This procedure does take into account hnite volume effects from three 
different sources: hrstly from the gluon propagator, secondly from the quark-gluon vertex 
via its ghost content in r 2 and thirdly effects generated by the quark-SDE itself. The 
parameters of the htted quark-gluon interaction are kept hxed, as are the renormalised 
current quark masses. In Fig. Elwe compare our results in the inhnite volume limit to the 
lattice data. We observe hnite volume ehects for momenta smaller than 1 GeV. The most 
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FIG. 6: The quark mass function at a given momentum Mm := M{p‘^ = 0.0924GeV^) plotted as a function of 
the current quark mass in the continuum and on tori with different volumes (staggered quarks in the left 

diagram, overlap quarks in the right diagram). 


pronounced effects occur in the infrared momentum region, where the continuum mass 
functions approach hnite values at vanishing momentum which are substantially larger 
than anticipated from the lattice data alone. We conclude that lattice data underestimate 
the quark mass that is dynamically generated by a substantial amount. The absolute size 
of this effect is approximately the same for all quark masses we have investigated, so the 
relative size becomes smaller for heavier quarks. Exactly as expected heavier quarks £t 
on hnite volumes most readily, but here we have deduced the size of the volume effects 
for all masses. 

Another interesting effect occurs for small momenta: the quark mass function in the 
chiral limit is actually larger than the one for the smallest nonzero current quark mass. 
This effect has also been observed directly from the quenched staggered lattice data . 
The effect is more pronounced for the staggered quarks, as can be seen by comparing the 
two diagrams of Fig. IHlat small quark masses. There we show the quark mass function 
at a given momentum. Mm := = 0.0924 GeV^), plotted as a function of the cur¬ 

rent quark mass. Apart from small quantitative differences, the qualitative behaviour of 
Min{m) is similar for both fermion formulations. On the compact manifold, one observes 
a (small) volume dependence for 2 < L < 3 fm, which becomes negligible for L > 3 fm. 
However, even then one is still far away from the inhnite volume results. This remaining 
difference stems from the Yang-Mills sector of the theory: the difference there between 
the torus and inhnite volume results (c/. Fig. |2)) has a direct impact on the ehective 
interaction in the quark-SDE and therefore on the results for the quark propagator. The 
size of this ehect depends on the continuum value of the exponent n (c/. ea. (ITT|l i. which is 
closely related to the quality of the approximation of the ghost-gluon-vertex in the ghost 
and gluon SDEs. There are indications that the ‘true’ exponent k may be closer to k ~ 0.5 
and therefore closer to the current lattice data than our value k, ~ 0.596 miillli. Thus 
we regard our result as an upper limit for the ehects that depend upon the volume. 
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L [fm] 

FIG. 7: The quark mass function for an up/down quark at a given momentum p = IGeV plotted as a function of 
the box length. 


C. A critical volume for chiral perturbation theory 

As discussed above, dynamical chiral symmetry breaking alone cannot occur on a finite 
volume. A nonvanishing current quark mass satisfying the condition 

Vm\{qq)\ » tt (23) 

has to be present as a ‘seed’ to trigger dynamical mass generation and the formation of a 
chiral condensate. Chiral perturbation theory, on the other hand, is built upon the chiral 
limit and therefore can only be applied on volumes large enough such that very small 
quark masses are feasible. The corresponding condition, discussed above, is 

L » ~ 1 fm. (24) 

^JTT 

We will now give an estimate of how large L has to be in practice. To this end we employ 
a current quark mass of the order of a typical up/down-quark mass, M{p^ = 2.9GeV^) = 
lOMeV , and determine the mass function M{p^) at = 1 GeV from solutions on tori with 
different volumes by linearly interpolating on the corresponding momentum grids. The 
result is plotted in Fig. 0 We clearly see that the quark mass function grows rapidly in 
the range 1.0 < L < 1.6 fm signalling the onset of dynamical chiral symmetry breaking. 
Above L = 1.6 fm, a plateau is reached. This picture does not change when we extract 
the mass function M(p^) at smaller momenta or when we employ even smaller quark 
masses. Thus a safe value for L should be at least 

Ly^PT ^ 1.6 fm. (25) 

This is a surprisingly small value in the light of the condition Eq. (El and so provides 
some justification for extending chiral perturbation theory to rather small volumes. 
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V. PROPERTIES OF QUARKS AND PIONS IN THE INFINITE VOLUME/CONTINUUM 

LIMIT 

A. The chiral condensate 

Having studied the volume dependence of the quark propagator in some detail, we now 
focus on and investigate the quark/meson sector employing our effective interaction 
which has been hxed by the lattice input. From the quark propagator S-^ in the chiral limit, 
we can determine the value of the (renormalisation point dependent) chiral condensate 
using 

-(??)(h^) := J , (26) 

where the trace is over Dirac indices. Its value is conveniently determined at a large renor- 
malisation scale, converted to the renormalisation point independent chiral condensate 
and then run down to /i = 2 GeV employing the quenched scale ^qcd — 0.225(21)MeV 
j^ . We then obtain the values 

-{m)TeriaM) = (252.6 ± b.OMeV)^, -{qq)Taaeredif^") = (253.0 ±5.0MeV)3, (27) 

which are in very good agreement with each other. The error given is an estimate of 
combined numerical and scale uncertainties. It is a quite interesting and satisfying result 
that the continuum chiral limit quark propagators of the staggered and overlap fermions 
agree very well with each other. It is furthermore interesting to compare our value for 
the chiral condensate with recent results using other methods on the lattice. Gimenez et 
al hnd the value (265 ± 27MeV)^ from an operator product expansion employing an 
0(a)-improved quenched Wilson action. Wennekers and Wittig quote (285±9MeV)^, 
determined from a quenched overlap action. Both values are in fair agreement with each 
other and with our result. McNeile recently obtained the value (259±27MeV)^ from 
a chiral Lagrangian with parameters hxed by lattice data employing Nf = 2 + 1 staggered 
sea quarks, thus indicating that unquenching effects in the chiral condensate may be small. 
This is in excellent agreement with the prediction in the SDE/BSE-approach |l7|, |21|. 


B. On the analytical properties of the quark propagator 


The analytic properties of the quark propagator can in part be read off from the corre¬ 
sponding Schwinger function 


a{t) = / d?x 


+7 


exp(!p ■ a:)cTsy (p^), 


(28) 


where asy are the scalar and the vector parts, resp ectiyely , of the dressed quark propa¬ 
gator. (This method has a long history, see e.g. 


resp ectiyely , oi tne aresse^ 

IE 13 11) IE iE ill and references 

therein). According to the Osterwalder-Schrader axioms of Euclidean held theory li 


this function has to be positive to allow for asymptotic states in the physical sector of the 
state space of QGD. Gonversely, positivity violations in the Schwinger function show that 
the corresponding asymptotic states (if present) belong to the unphysical part of the state 
space. Thus positivity violations constitute a sufficient condition for conhnement. Our 
results for the Schwinger function of the chiral limit quark propagator in the staggered 
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FIG. 8: The logarithm of the Schwinger function ln(|(j(f)|) of the chiral limit quark propagator as a function of 
time. Shown are results for the staggered and overlap quark, (i) employing the ansatz Eq. for the quark-gluon 
vertex and (ii) substituting the Ball-Chiu construction, Eq. (HU, for The respective curves for the staggered 
and overlap quarks are indistinguishable in the plot. We compare the results to fits of the function Eq. iHU. 


and overlap formalisms are shown in Fig. |H1 The Schwinger functions of the two formu¬ 
lations agree extremely well and are indistinguishable in the plot. Let us hrst discuss the 
result obtained with our simple vertex construction involving only the vector part 7 ^ of 
the vertex. The cusp at t = 6.76 GeV“^ = 1.33 fm indicates a node in the Schwinger 
function corresponding to positivity violations at a scale in rough agreement with the size 
of hadrons. An excellent £t to the Schwinger function is obtained using the form 

a{t) = \boexp{-bit) cos{b 2 t + h) \ , (29) 


which corresponds to a pair of complex conjugate poles of the propagator in the timelike 
momentum plane. These poles correspond to a ‘quark mass’ given by m = 61 ±ib 2 , which 
in our case is m = 516(20) ± i 428(20) MeV. Taken at face value this means that the 
lattice quark propagator is confined (this is also the conclusion drawn by Bhagwat et 
al. 0]). However, there is a caveat: it has been shown in Ref. that the presence of 
a sufficiently strong scalar part in the quark-gluon vertex can have a strong influence on 
the analytic structure of the solution of the quark-SDE. Indeed if we employ the Ball- 
Chiu construction, Ea. () 22 j) . instead of 7 ^ we obtain an exponentially decaying Schwinger 
function denoted by the dash-dotted line in Fig. |S1 Such a function corresponds to a 
positive definite quark propagator with a pole on the real axis at m = 632(20) ± i0{2) 
MeV (within numerical accuracy). This shows that the analytic structure of the quark 
propagator depends strongly on the details of the structure of the quark-gluon vertex and 
one cannot make definite statements from fitted interactions alone. First attempts on the 
lattice ^ as well as in the SDE/BSE approach H H US have been made to study 
the tensor structure of the quark-gluon vertex in more detail. Still, more effort is needed 
before definitive conclusions about the analytic structure of the quark propagator are in 
sight. 


C. Pion mass and decay constant from the lattice interaction 

Finally we have determined the mass and the decay constant of the pion employing the 
effective quark-gluon interaction fitted to the lattice data. The pion is described by the 
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P(MeV) 

staggered 

4.1 ±0.3 

138.4 

83.5 

overlap 

4.1 ±0.3 

138.7 

84.3 

experiment 


138.5 

92.4 


TABLE III: Results for the renormalised current up/down quark mass, the mass of the pion and the pion decay 
constant employing the quark-gluon interaction fitted to the lattice data. The evolution of the quark mass has 
been performed using ^qcd = 0.225(21)MeV 1^ . 


homogeneous Bethe-Salpeter equation (BSE) 

/ d'^k 

7^ fc; P) [Sik^W’ik P)S{k_)]^, (30) 

where 


r-(p; P) = 75 P) - P) - P) - H-{p; P )], (31) 


is the Bethe-Salpeter amplitude of the pion, K is the Bethe-Salpeter kernel and the 
momenta = k + and k^ = k + — 1)P are such that the total momentum P = 

— /c_. All physical results are independent of the momentum partitioning parameter 
^ = [0,1]. The crucial link between the meson bound states and their quark and gluon 
constituents is provided by the axial vector Ward-Takahashi identity. It relates the quark 
self energy to the quark-quark interaction kernel in the BSE and thereby guarantees the 
Goldstone nature of the pions and kaons . In our case the kernel is given by 


Ko,f3-,5^{p,k-,P) [ 7 ^]^^ [-f^]g^tf,^{k)Z. 


Z{k) T,{k) V2{k) V^{k) 


(32) 


where is a transverse projector in momentum space and the flavour content of the 
kernel has been suppressed. To determine the pion mass and wave functions we explicitly 
solve the quark-SDE in the complex momentum plane, thereby providing the necessary 
input for the BSE. The BSE is then solved as an eigenvalue equation for eigenvalue one, 
which provides the wave functions and the pole mass of the pion (all technical details 
of such a calculation are discussed in detail in Ref. ^3). From the (normalised) wave 
function the pion decay constant is hxed by 


Q r ph 

A = — Tr, j ^r-(fc, -P)S{k + P/2)75/(fc - P/2). (33) 

Our results for the pion mass, the corresponding renormalised current quark mass in the 
MS-scheme and the pion decay constant are given in Table ITTTl The current quark mass 
has been determined from the quark-SDE at a large renormalisation point, converted into 
the MS'-scheme and sub seq uently evolved to p = 2 GeV employing the same scale as for 
the chiral condensate, = 0.225(21)MeV [^. The errors given in Table imi are an 

estimate of numerical and scale uncertainties. The resulting current quark mass is in the 
ballpark of the values quoted by the particle data group |53| . 
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Probably the most interesting observable is the pion decay constant, which directly 
reflects the deflciency of the quenched lattice calculation compared with the real world. 
Both lattice formulations underestimate the experimental value by rou ghly ten percent. 
This margin is much smaller than the thirty percent estimated in Ref. where finite 
volume effects were not taken into account. 


VI. SUMMARY 


In summary, we have investigated the properties of a quenched lattice-QCD quark 
propagator from staggered quarks [l^ and from an overlap quark action Q in Landau 
gauge. We employed a coupled set of Schwinger-Dyson equations for the ghost, gluon and 
quark propagators on compact manifolds and in the infinite volume/continuum limit to 
study finite volume effects in the propagators and to determine some aspects of dynamical 
chiral symmetry breaking on a torus. We constructed a model for the quark-gluon vertex 
such that two sets of staggered and overlap lattice quark propagators are reproduced by 
the quark gap equation on their respective manifolds. 

Comparing results on different volumes with the infinite volume/continuum limit, we 
found sizeable quantitative but not qualitative differences at small momenta. The contin¬ 
uum mass functions approach finite values at vanishing momentum which are substantially 
larger than anticipated from the lattice data alone. Thus lattice simulations may underes¬ 
timate the amount of dynamical quark mass generation in the infrared by as much as 100 
MeV. The absolute size of this effect is approximately the same for all quark masses we 
have investigated, so the relative size becomes smaller for heavier quarks. As a by-product 
of this investigation we observed that the bare quark masses in the staggered and overlap 
formulations are related by a simple multiplicative factor. No additional additive correc¬ 
tions were necessar y, w hich is a signal that taste symmetry violations in the staggered 
action used in Ref. l3| are negligible 

We also assessed the effects of small volumes on dynamical chiral symmetry breaking. 
Employing a fixed small current quark mass we decreased the volume of the box until we 
found clear signals of chiral symmetry restoration. These signals occur at the surprisingly 
small box length of 


L ~ 1.6 fm, (34) 

which constitutes a minimal box size below which chiral perturbation theory cannot safely 
be applied. 

With the quark-gluon interaction fixed by the lattice data, we then determined the 
properties of the quark propagator and pions in the infinite volume/continuum limit. We 
found a chiral condensate of 


\{qq)\^^^ = {253 ±5 MeVf (35) 

which compares favourably with values determined with other methods on the lattice. 
Unfortunately nothing can be said about the analytic structure of the quark propagator. 
We showed, that the method of fitting a quark-gluon interaction to lattice propagator 
data is not sufficient to pin down the relative strength of the various tensor components 
of the vertex, which in turn are necessary to derive reliable statements on the analytic 
structure of the quark propagator. Finally we have determined the pion mass and decay 
constant employing a rainbow-ladder truncation of the Bethe-Salpeter equation, which 
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incorporates the fitted interaction. We obtained the renormalised up/down quark masses 

= 4.1 ± 0.3 MeV, (36) 

which are in the ballpark of the values given in the particle data book |^. The pion decay 
constant is roughly ten percent smaller than the experimental value. This indicates that 
quenching effects in the light meson sector are not too large in agreement with previous 
hndings |l7| . 

We have shown how the SDE/BSE approach, once matched to lattice results on hnite 
volumes with appropriate manifolds, can reliably determine inhnite volume/continuum 
predictions for all quark masses, including those of the real world close to the chiral limit; 
a limit not directly accessible on the lattice. 
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